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SUMMARY 

Several computational levels of progressive sophlstlcatlon/slmpllf Icatlon 
are described to computationally simulate composite sandwich hygral, thermal, 
and structural behavior. 

The several computational levels of sophistication Include: (1) three- 

dimensional detailed finite element modeling of the honeycomb, the adhesive 
and the composite faces; (2) three-dimensional finite element modeling of the 
honeycomb assumed to be an equivalent continuous, homogeneous medium, the 
adhesive and the composite faces; (3) laminate theory simulation where the 
honeycomb (metal or composite) Is assumed to consist of plies with equivalent 
properties; and (4;) derivations of approximate, simplified equations for ther- 
mal and mechanical properties by simulating the honeycomb as an equivalent 
homogeneous medium. The approximate equations are combined with composite 
hygrothermomechanlcal and laminate theories to provide a simple and effective 
computational procedure for simulating the thermomechanlcal/thermostructural 
behavior of fiber composite sandwich structures. 


INTRODUCTION 

The use of fiber composites In space applications Is Increasing In a 
variety of structural configurations. Sandwich structural configurations pro- 
vide an effective application of fiber composites. The faces.of the sandwich 
resist loads by membrane action which Is one of the most structural ly-efflclent 
use of fiber composite thin laminates.. In addition,, composite sandwich struc- 
tures can be designed to meet very close thermal distortion tolerances such as 
those required for communication satellite antennas and reflectors. 
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Designs to meet close thermal distortion tolerances are developed by 
detailed thermal and structural analyses. These analyses require thermal and 
mechanical properties of the composite sandwich as well as temperature and 
moisture effects on these properties. The thermal properties required Include 
heat capacity and thermal conductivities and thermal expansion coefficients In 
the plane and through the thickness of the sandwich. The corresponding mechan- 
ical properties Include normal and shear moduli, and Poisson's ratios. Thermo- 
mechanical and corresponding thermostructural properties of sandwich components 
can, In principle, be measured experimentally. This Is practical only for one 
or at most a few sandwich configurations, limited temperature/moisture condi- 
tions and generally for selected properties. 

The alternative Is to computationally simulate the thermomechanical 
behavior of composite sandwich structures so that all the properties required 
for thermal and mechanical/structural analyses can be predicted. Recent 
analytical studies at Lewis Research Center focused on developing computational 
methods for simulating the thermomechanical behavior of composite sandwich 
structures. These methods use analyses with several levels of progressive 
sophlstlcatlon/slmpl If Icatlon In conjunction with composite hygrothermomechan- 
Ical theory. The objective of this paper Is to describe these computational 
simulation methods and summarize results obtained therefrom. 

The several computational levels of sophistication Include: (1) three- 

dimensional finite element modeling of the honeycomb, the adhesive and the 
composite faces; (2) three-dimensional finite element modeling of the honeycomb 
assuming an equivalent homogeneous medium, the adhesive and the composite 
faces; (3) laminate theory simulation where the honeycomb (metal or composite) 
Is assumed to be of plies with equivalent properties; and (4) approximate, 
simplified equations for simulating the honeycomb thermal and mechanical prop- 
erties with an equivalent homogeneous medium. These levels of sophistication/ 
simplification have been packaged Into a procedure which Is embedded In a com- 
posite mechanics computer code (ref. 1) streamlined for the computational simu- 
lation of composite sandwich hygral, thermal, and structural behavior. The 
steps for developing such a procedure and Its subsequent embedment. In a com- 
puter code are described In outline form. 


THREE-DIMENSIONAL DETAILED FINITE ELEMENT MODEL 

The three-dimensional detailed finite element simulation Includes simula- 
tion of the (l) core, (2) adhesive, and (3) faces. The core Is simulated using 
three plate elements through the core thickness In order to determine equiva- 
lent mechanical and thermal properties and possible core buckling. The finite 
element model of a core generated by MSC/NASTRAN Is shown In figure 1. The 
core for this model Is 0.375 In. high. The model consisted of 380 nodes, 375 
elements, and 2280° of freedom (DOF). The model was used to determine nine 
equivalent mechanical properties (^exx* ^eyy* Eg ZZ , ^exy* ^eyz* ^ezx f ^exy 
«eyz* v ezx) and six equivalent thermal properties (c*exx» a ey y, <* ezz , K exx. 

K ey y , K ezz ). The notation Is as follows: E denotes normal (Young's) modulus, 

G denotes shear modulus, v denotes Poisson's ratio, a denotes thermal 
expansion coefficient, and K denotes thermal heat conductivity. The sub- 
script e denotes core equivalent property while the subscripts x, y, and 
z denote core coordinate reference plane (first) and direction (second). The 
material properties used with the finite element model are summarized In table I. 
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The equivalent mechanical properties of the core are determined by loading 
the model with Imposed displacements In each direction at one-end (or face) 
while fixing the opposite end (or face). The desired property Is determined 
by using corresponding mechanics of materials equations. For example, E e22 
= N CX x^x/^x u where N C xx the reaction force due to Imposed displacement 
u at. x = lx, lx Is the length of the core In the x direction. Ax Is the 
core projected area on a plane with normal along the x direction and u Is 
the Imposed displacement. The Poisson's ratio v CX y = 2vl x /uly where v Is 
the average of all nodal displacements at the y = 6 and y = ly faces and u 
Is the Imposed displacement. The remaining mechanical properties are determined 
In a similar manner. The thermal properties are also determined In a similar 
manner, except that the Imposed displacements are replaced with appropriate 
thermal conditions. 

Deformed finite element plots superimposed over the undeformed models are 
shown in figures 2 to 7 for the nine different mechanical properties. These 
figures are Included to show that the honeycomb wall did not undergo large 
deformations In any of the different loading conditions so that no linear 
solution Is required. Displacements were Imposed In two planes, for example, 
xy and yx (figs. 5 to 7) In order to determine the equivalent shear moduli 
and their respective symmetries. The equivalent mechanical properties of the 
core determined from these cases are summarized In table II for two cell sizes 
(1/4 and 1/8 In.). As can be seen: (1) the normal moduli are small because 

the honeycomb deforms like an accordion In these directions, and (2) some 
shear moduli are unequal and Poisson's ratios do not satisfy the well known 
reciprocity rules. That Is, the y property Is not equal to the x property, 
for example. This Implies that the behavior of the core Is load-direction 
dependent when the mechanical properties are determined from the simple, basic 
mechanics of materials definitions. It Is Important, therefore. In measuring 
and reporting these properties to Include the specific plane and direction In 
order to use them properly In simulated core or structural sandwich analysis.. 
The unequal shear moduli and Poisson's ratios demonstrate that the honeycomb 
core behaves like an equivalent homogeneous anisotropic solid. 

The computational simulation of the entire sandwich Is parallel to that 
of the fabrication procedure. First the core Is simulated as already described. 
Second the adhesive layers are simulated with a single layer of solid finite 
elements on both sides of the core. Third, the composite faces are simulated 
with a single layer of solid elements over the adhesive on each side. An 
xploded schematic of the MSC/NASTRAN finite element model Is shown In figure 8. 

This model consists of 760 nodes, 711 elements, and 2850 DOF. The material 
properties used are those In table I. 

The equivalent properties of the sandwich are determined by loading the 
model with Imposed displacements as was done for the core. Typical results 
obtained are summarized In table III under , detal led model. It can be seen 
from these results and referring to table I that the In-plane properties of 
the sandwich are mainly controlled by the faces while the flat-wise (through- 
the-thlckness) properties are controlled by the core. 


THREE-DIMENSIONAL FINITE ELEMENT MODEL WITH EQUIVALENT CORE 

This approach Is the first step In the progressive simplification for 
simulating the composite honeycomb sandwich using laminate theory. The finite 
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element model for this simulation consisted of four layers of solid elements 
for the core, one single layer of solid elements for the adhesive on each side 
of the core, and two layers of solid elements for each face. A computer 
schematic of the finite element model used for determining the structural 
properties Is shown In figure 9. The model consists of 2475 nodes, 1960 ele- 
ments, and 7425 OOF. The material properties for the model were those In 
table I for the adhesive and the faces. Those for the honeycomb core were the 
equivalent properties determined from the three-dimensional detailed finite 
element model listed In table II. The properties of the grouping of plies 
represented by the solid elements for the faces were determined using a com- 
posite mechanics code ICAN (ref. 1). Each grouping consisted of two (+60, 

-60, 0) layers to simulate the six ply [±60, O 2 , +60] face sheet In order to 
have a minimum thickness laminate with In-plane Isotropic behavior for mechan- 
ical and thermal properties. A 40-1n. -square sandwich with the same thickness 
and same modeling detail was used for determining the thermal expansion coeff- 
icients In order to obtain uniform thermal expansion In all planes. This fin- 
ite element model was loaded to determine the sandwich structural and thermal 
properties as was done for both the honeycomb core and the three-dimensional 
detailed finite element model. Results obtained from this simulation are 
described later In the results and comparison section. 


LAMINATE THEORY 

Laminate theory Is the second level In the progressive simplification for 
simulating structural sandwich behavior. In the studies summarized herein the. 
sandwich was simulated using laminate theory as follows: (1) four plies, for 

the core, (2) one ply for each of the adhesive layers, and (3) six plies for 
each of the faces. The desired structural and thermal properties are directly 
determined from the analysis and are routine output properties from the ICAN 
code. 


The input properties for this simulation are those used for composite 
micromechanics and most of them are available In the ICAN resident data bank. 
The equivalent properties for the honeycomb were modeled using assumed fiber/ 
matrix properties to match those determined from the core three-dimensional 
detailed finite element model. It Is Important to note that In this simulation 
the number of plies used for the faces Is equal to that In the actual sandwich. 
The number of plies for the adhesive and the core Is the choice of the user. 
These can be selected to account for nonuniform temperature and moisture pro- 
files. The hygrothermal effects on the sandwich structural and thermal prop- 
erties can readily be predicted using laminate theory which accounts for these 
effects. Properties obtained from this simulation using ICAN are summarized 
In the next section. 


RESULTS, COMPARISONS AND DISCUSSION 

The results obtained from the three different simulation methods are sum- 
marized In table III. As can be seen the results from all three simulation 
methods are generally In close agreement except for some Poisson's ratios as 
explained later. Considering the three-dimensional detailed finite element 
model as being the most accurate. It can be seen that laminate theory can 
reasonably accurately predict the structural and thermal properties of sandwich 
structures. This Is a significant finding since composite mechanics computer 



codes can be used to simulate structural and thermal sandwich behavior and 
conduct parametric studies. Including hygrothermal effects, In a near-routine 
and very cost-effective manner. Typical thermal expansion coefficient results 
obtained In this manner (using ICAN) are plotted In figures 10 and 11 for com- 
posite sandwiches with different laminate configurations and with different 
honeycomb material suitable for applications to communication satellites. The 
near-zero In plane thermal expansion coefficients as well as the hygrothermal 
effects on these coefficients are Illustrated. The apparent discrepancies are 
mainly due to the effects of the free edges In the relatively small size model 
for the three-dimensional detail finite element. These same effects also 
Influence the a zz thermal expansion coefficient. The model-size free-edge 
effects were evaluated using the three-dimensional finite element model with 
homogeneous core. A finite element model size of 40 by 40 In. was necessary 
to practically eliminate the free-edge effects. 


APPROXIMATE EQUATIONS FOR EQUIVALENT CORE PROPERTIES 

A set of simple equations are presented for predicting homogeneous-core 
equivalent honeycomb properties. These equations are the third approach to 
the progressive simplification of composite sandwich panels. The equations 
can be used to predict approximate properties for a honeycomb core made from 
any material , for example, metal, nomex, or composite. The equations are 
derived with the aid of figure 12. Note the figure reference coordinate axes 
and the schematics showing the stresses In the honeycomb core wall. The sub- 
scripts used In the notation for the core properties refer to this figure. 

The equations were derived assuming that the honeycomb walls resist/react 
stress through membrane (Inplane) action only. The rationale for assuming 
only Inplane action Is that the sandwich faces will prevent excessive bending 
and/or linkage-type behavior of the honeycomb walls. 

The resulting equations are summarized In table IV together with their 
respective predictions. Included are also predictions from finite element 
simulations for (1) a single cell with faces to prevent core wall bending and 
without faces, and (2) the three-dimensional detailed finite element model for 
the core only. The coefficients In the approximate equations were slightly 
adjusted from those derived for the mechanistic models so that predictions 
correlate with the finite element simulations. It Is worth noting that the 
approximate equations depend only on core wall thickness (t), cell size (w) 
and core wall material properties. It Is also worth noting that core proper- 
ties E exx and E eyy are practically "zero" and u exy Is approximately 
unity for the core when It Is not restrained. The stress-strain states used 
to determine these unrestrained core properties require the core to respond as 
an Integrated linkage system because of the large bending displacements of the 
core wall which are very thin (0.0007 In.). 

A procedure was developed to use these approximate equations In conjunc- 
tion with composite mechanics In order to simulate sandwich thermal/structural 
behavior. The steps of this procedure are as follows: (1) obtain the core 

material properties from suppliers, (2) use the approximate equations to 
determine equivalent homogeneous core properties, (3) select fiber/matrix 
properties with the aid of composite micromechanics equations (ref. 2) to 
reproduce the core equivalent properties, (4) use the micromechanics equations 
to predict all ply properties needed for laminate analysis, (5) select the 
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number of plies to represent the adhesive and the core, and (6) use laminate 
theory (refs. 3 to 5) to predict the desired properties. 

The procedure just outlined Is embedded In the ICAN computer code. In 
addition, ICAN was appropriately modified and packaged as a stand-alone port- 
able code for the computational simulation of structural sandwich hygrothermo- 
mechanlcal behavior. The resulting code Is Identified as ICAN/SCS for Inte- 
grated Composite Analyzer for Structural-Sandwich Computational Simulation. 
ICAN/SCS output Includes material cards for solid and plate finite element 
structural analysis as well as ply and Interply stress analysis. The ply and 
interply stress analysis of the core has to be carefully Interpreted since the 
honeycomb is replaced with an equivalent homogeneous, continuous anisotropic 
solid. 


SUMMARY 

The thermal and structural behavior of composite sandwich panels with a 
honeycomb core has been computationally simulated using several computational 
levels of progressive sophisticatlon/slmplif ication. These levels include: 

(1) three-dimensional detailed finite element modeling, (2) three-dimensional 
finite element modeling assuming a homogeneous core, (3) laminate theory, and 
(4) simple equations for predicting the equivalent properties of the honeycomb 
core. A procedure was developed and embedded In a composite mechanics computer 
code which can be used to conduct parametric studies In order to determine 
"optimum" composite sandwich configurations for specific applications. The 
procedure developed makes it possible and computationally effective to evaluate 
composite sandwich behavior at the global, local, laminate, ply, and micro- 
mechanics levels when the composite sandwich is subjected to hygral, thermal 
and mechanical loading environments. The cross correlation at the various 
computational levels of progressive simplification provide credence that the 
sandwich behavior so simulated Is representative and consistent with the 
assumed physics. 
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TABLE I. - PROPERTIES USED AS INPUT TO THE SIMULATION MODELS 


Property 

Units 

Symbol 

Direct. 

Aluminum 

core 

Unidirection 

composite 

(Ply 3 ) 

(T-300/IMHS) 

Composite 

face 

sheet 

Adhesive 
1 ayer 

Normal 

Mpsi 

£ 

XX 

10.0 

19.4 

7.312 

0.3 

modulus 



yy 

10.0 

1.193 

7.312 

.3 




ZZ 

10.0 

1.193 

1.407 

.3 

Shear 

Mpsi 

G 

xy 

3.85 

0.552 

2.783 

.1034 

modulus 



yz 

3.85 

.331 

0.442 

.1034 




ZX 

3.85 

.552 

.442 

.1034 

Poisson 1 s 


\> 

xy 

0.312 

.26 

.3135 

.45 

ratio 



yz 

.312 

.424 

.3135 

.45 




ZX 

.312 

.26 

.275 

.45 

Thermal 

uin./in./’F 

a . 

XX 

12.6 

- .1732 

.9208 

57.02 

expansion 



.yy 

12.6 

15.26 

.9208 

57.02 

coefficient 



ZZ 

12.6 

15.26 

21.04 

57.02 

Thermal heat 

Btu-in./hr 

K 

XX 

6.166 

2.420 

1.225 

.0087 

conductivity 

in.^/'F 


yy 

6.166 

.0297 

1.225 

.0087 




zz - 

6.166 

.0297 

.0297 

.0087 

Density 

lb/ in. 3 

P 

— 

.097 

.056 

.056 

.04 

Heat capacity 

Btu/lb 

C 

— 

.22 

. .1951 

.1951 

.25 


a The x-x direction is taken parallel to the fiber direction, y-y perpendicular to 
it, and z-z through the thickness. 


TABLE II. - CORE EQUIVALENT MECHANICAL PROPERTIES 
DETERMINED BY USING THREE-DIMENSIONAL DETAILED 
FINITE ELEMENT MODEL 


Property 

Units 



Aluminum honeycomb 
cell size 

1/4 in. 

1/8 in. 

Normal 


E 

XX 

0.00178 

0.01423 

modulus 



yy 

.00175 

.01401 


Mpsi 


ZZ 

.07541 

.15083 

Shear 

Mpsi 

G 

xy 

.00797 

.01560 

modulus 

psi 


yx 

.00125 

.01002 


Mpsi 


yz 

.01065 

.02147 




zy 

.01065 

.02147 


Mpsi 


ZX 

.01561 

.03185 




xz 

.01775 

.03572 

Poisson's 


V 

xy 

1.029 

1.029 

ratio 



yx 

1.024 

1.024 




yz 

0 

0 




zy 

.311741 

.311692 




ZX 

.311633 

.311680 




xz 

0 

0 

Volume ratio 




.007467 

.014934 



























TABLE III. - SUMMARY OF RESULTS FROM THREE DIFFERENT 
SIMULATION LEVELS OF PROGRESSIVE SIMPLIFICATION 


Property 

Units 

Symbol 

Direct. 

Three-dimensional finite 
element 

Laminate 

theory 

Detailed 

model 

Homogeneous 
core model 

Normal 

Mpsi 

E 

XX 

1.0 

1.0 

1.0 

modulus 



yy 

1.0 

1.0 

1.0 




ZZ 

.09 

.09 

.09 

Shear 

Mpsi 

G 

xy 

.39 

.38 

■ .38 

modulus 



y * 

— 






yz 

.071 

.071 

.071 




zy 







zx 

.075 

.074 

.076 




xz 




Poisson's 


V 

xy 

.396 

.315 

.314 

ratio 



yx 

.348 

.315 

.314 




yz 

.121 

.043 

a 0 




zy 

.073 

.005 

a 0 




ZX 

.030 

.004 

a 0 




XZ 

.118 

.045 

a 0 

Thermal 

uin./in./'F 

a 

XX 

1.312 

1.21 

1.21 

expansion 



yy 

1.378 

1.21 

1.21 

coefficient 



ZZ 

17.791 

14.6 

13.6 

Thermal heat 

Btu-in./hr 

K 

XX 

.191 

.191 

.190 

conductivity 

in.2/°F 


yy 

.181 

.183 

.182 




ZZ 

.043 

.043 

.041 


a Were set to "zero" in the laminate theory. 






















TABLE IV. - SUMMARY OF EQUIVALENT CORE PROPERTIES OF 1/4 IN. ALUMINUM HONEYCOMB 


Property 

Units 

Symbol 

Modulus 

ksi 

Eexx 



Eeyy 

^ezz 

p ex y 

b eyz 

b ezx 

Poisson's 

ratio 


v exz 

v eyz 

y ezx 

Thermal 

expansion 

coefficient 

uin./in./°F 

“exx 
“e yy 
“ezz 

Thermal heat 
conductivity 

Btu/in./hr 

in. 3 /°F 

’ K exx 

£eyy 

K ezz 

Heat capac- i 
ity density 

Btu/lb 

c e 

Oensity 

lb/in. 3 

p e 


Approximate equations 

Finite element analysis 

Faces 

restrained 

Faces 

free 

Faces 

restrained 

Faces 

free 

Single cel 1 model 

Three-dimensional 
detail model 

Faces 

restrained 

Faces 

free 

1.0 (t/w)E hxx 

1.0 (t 3 /4w)E hxx 

29.6 

0.003 

28.5 

0.002 

0.0 


1.5 (t 3 /4w)E hxx 

44.4 

.003 

43.2 

.002 

.0 

3 (t/w)E h * 

3.0 (t/w)E hzz 

89.1 

89.1 

69.7 

69.0 

75.4 

(3/4) (t/w)G hxy 

(t/w) G h “ 

22.3 

7.4 

23.9 

6.9 

7.8 

1 .0(t/w)G hxz y 

1.0 {tlvi)Gfr 77 

11.4 

11.4 

10.5 

8.4 

10.6 

1.5(t/w)G hzx 

1.0 (t/w)G hzx 

17.1 

11.4 

14.5 

8.0 

7.8 

y hxz 

27 yj 

1.15 

.34 

.34 

1.0 

1.03 

v h 7 V 

0 

.30 

0 

.21 

.00 

.0 

v hzx 

v hzx 

.30 

.30 

.25 

.30 

\ 

.31 

“hxx 

°hxx 

12.6 

12.6 

12.6 

12.6 

12.6 

“hxx 

°hxx 

12.6 

12.6 

12.6 

12.6 

12.6 

“hzz 

“hzz 

12.6 

12.6 


12.6 

12.6 

(4/3) (t/w)K hxx 

(4/3){t/w)K hxx 

0.023 

.023 

— 


.027 

1.5(t/w)K hxx 

1.5(t/w)K hxx 

.026 

.026 


— 

.017 

3(t/w)K hzz 

3.0(t/w)K hzz 

.052 

.052 



.046 

3(t/w)C h 

3.0(t/w)C k 

.002 

.002 

— 


— 

3 ( t/w) oh 

3.0(t/w)oh 

.00084 

.0008 

— 


.0007 






























X 

Figure 3. - Deformed shape due to imposed Y- displacement for equivalent 
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(b) X on Y- face. 


Figure 5. - Deformed shape due to imposed displacement 
for equivalent G exy and G eyx . 




Figure 6. - Deformed shape due to imposed displacement for 
equivalent G exz and G ezx . 
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Figure 7. - Deformed shape due to imposed displacement 
for equivalent G eyz and G e2y . 



Figure 8. - 3-D detailed firfite element model of composite sandwich with 
honeycomb core. (All sandwich components are modeled with solid 
finite elements except honeycomb core. ) 





Figure 9. - 3-D finite element model of composite sandwich with equivalent homogeneous core. 
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(c) CTEXX-CTEYY versus temperature 
for aluminum honeycomb core (1/4- 
in. cells) with P-75/R934(0/±60)S 
composite faces. 


(d) CTEZZ versus temperature for alum- 
inum honeycomb core (1/4-in. cells) 
with P-75/R934(0/±60)S composite 
faces. 


Figure 10. - Thermal expansion coefficients of composite sandwich panels with graphite 
fiber/epoxy faces and aluminum honeycomb core. 
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Figure 12. - Schematic for deriving approximate equations for equivalent honeycomb properties. 
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